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Abstract. The spatio-temporal dynamics of the deformation of a vibrated plate is measured by a high 
speed Fourier transform profilometry technique. The space-time Fourier spectrum is analyzed. It displays 
a behavior consistent with the premises of the Weak Turbulence theory. A isotropic continuous spectrum 
of waves is excited with a non linear dispersion relation slightly shifted from the linear dispersion relation. 
The spectral width of the dispersion relation is also measured. The non linearity of this system is weak as 
expected from the theory. Finite size effects are discussed. Despite a qualitative agreement with the theory, 
a quantitative mismatch is observed which origin may be due to the dissipation that ultimately absorbs 
the energy flux of the Kolmogorov-Zakharov casade. 

PACS. XX.XX.XX No PACS code given 



The kinetic Weak Turbulence Theory (WTT) describes 
the long time statistics of a large number of waves ex- 
changing energy through a weak non linear coupling. It 
has been applied to a vast number of systems such as 
Alfven waves in solar wind, oceanic surface or internal 
waves, optics in non linear media, Bose-Einstein conden- 
sates... (see for a collection of examples). It predicts 
in particular a Kolmogorov-Zakharov energy cascade as 
a stationary solution of the non equilibrium case. Energy 
is transferred conservatively from the large forcing scales 
down to small scales (at which it is dissipated) through 
multiple waves resonances. In some cases, an inverse cas- 
cade can be observed if the wave action is conserved as 
well. This phenomenology shares many similarity with hy- 
drodynamical turbulence (either in 2D or 3D) so that in- 
formation provided by the WTT could be of interest in 
the quest of the still missing statistical theory of fluid tur- 
bulence. 

The main hypotheses of the WTT are (i) a infinite 
size system and (ii) weak non linearity of the waves, (i) 
ensures that frequencies are not discrete so that resonance 
conditions can be fulfilled, (ii) is required in the theory so 
that to have scale separation between linear (fast) time 
scales and non linear (slow) time scales. In laboratory ex- 
periments, (i) can obviously not be achieved rigorously 
and waves have finite amplitudes so that the validity of 
both (i) and (ii) can be discussed. This is especially the 
case for gravity surface waves in which both hypotheses 
could be non valid simultaneously [3lll]. Another concern 
in experiments is dissipation. The theory requires an ad- 
ditional hypothesis: (iii) the energy sinks are supposed to 
act at scales widely separated from the forcing so that en- 
ergy flows conservatively in a inert ial range of scales. This 
condition is often difficult to achieve in experiments, let 



alone in numerics. This may also alter the statistics of the 
waves [5 . Additional hypotheses are (iv) locality of inter- 
actions in /c-space and (v) phases of interacting modes are 
random. 



In many cases, measurements are difficult to imple- 
ment as in measurements of solar winds or in confined 
plasmas. So far measurements have been restricted mostly 
to spectral exponents. The measurements are often per- 
formed at a single point in space so that only temporal 
Fourier spectra can be accessed. Here we use a laboratory 
experiment which setup is sufficiently simple so that both 
the temporal and the spatial evolution of wave turbulence 
can be measured: wave turbulence in a thin elastic plate. 
This system has interests in its own (cymbals, vibration of 
structures in aeronautics) but also as a toy model of wave 
turbulence. The WTT has recently been applied to this 
system [6 and predicts 4- wave resonances. Experiments 
show a disagreement between the measured single point 
spectrum and the theoretical prediction [71[8]. The appli- 
cation of a high speed profilometry technique provides a 
measurement of the dynamics of the waves resolved both 
in time and space over a significant portion of the plate [9j . 
It is possible for the first time to get a measurement of the 
spatio-temporal structure of wave turbulence and thus to 
test the WTT on more elaborate predictions as just single 
point spectral exponents. 



In the following, we recall previous theoretical and ex- 
perimental results, then we present shortly the experimen- 
tal setup. We analyze the statistics of the wave amplitude 
and the various Fourier spectra of the deformation. Finally 
we discuss the issue of finite size effects. 
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1 Summary of previous results 



1.2 Experimental results 



1.1 Theoretical results 

The evolution equations for the deformation of a thin elas- 
tic plate are the Foppl-Von Karman equations: 
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where is the normal displacement and E 2 x 10^^ Pa 
is the Young modulus, a :^ 0.3 is the Poisson ratio, and 
p :^ 8000 kg.m~^ is the specific mass (values are typical 
for stainless steel). A is the 2D Laplacian, and {a, 6} = 

b. X is the stress function 
and is related to the stretching of the plate [lOJ. 

These equations support linear dispersive waves (Ihs 
terms in ([T])) whose linear dispersion relation (LDR) is 
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For large amplitudes, the wave equation exhibits cubic 
nonlinearities (rhs of ([T]) and (j2j)). These nonlinear equa- 
tions are believed to lead to wave turbulence [6 . The 
WTT has been applied to the case of thin elastic plates 
by Diiring et al. [6 . Starting from ([T]) and (|2]), the waves 
are supposed to be weakly non linear. This enables a 
slow enough transfer of energy among waves so that the 
wave amplitudes evolve on much longer time scales than 
the linear period ([3]). In that case a stationary out of 
equilibrium regime can be exhibited in which energy is 
transferred through 4-wave resonances from large scales 
to small scales. In the following, we refer to "space spec- 
trum" for a power spectrum density £^(k) using the k (or 
k = ||k||) variable, to "time spectrum" for a power spec- 
trum density E{uj) using the uj variable (or / = uj/27t) 
and "space time spectrum" for £^(k, cj) (or E(k^uj)). We 
add C, oi V {v = dtC) subscript to E{. . .) to specify 
which quantity is used (deformation or normal velocity). 
The WTT predicts the following space spectrum for the 
displacement 



E^{k) = C 



ln^/^(feVfe) 
(1-^2)1/6 ^/E/^ 
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where P is the energy fiux per unit mass (here we use the 
definition of Connaughton et al. [11 ), C a dimensionless 
number, and k'^ is a cutoff wave number. Using (j3j) to 
change variables from k to the theoretical prediction 
for the velocity time spectrum is then 



(1 - £72)2/3 
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where oo^ = uoLik'^) and C is a number. 



The motion of thin steel plates has been studied by Boudaoud 
et al. [7 and Mordant [8 using single or differential mea- 
surement of the velocity with a laser vibrometer. Single 
point measurements provide only the time spectrum. The 
experimental velocity spectra have be shown not to obey 
the analytical prediction (|5|) neither in the frequency scal- 
ing nor in the energy fiux scaling. The observed scaling 
is Ey{uj) - pO.7^-0.6 ^^^Yier than Ey{uj) - P^/^cjO. A 
continuous spectrum was observed at the highest forcing 
and differential measurements are consistent with a su- 
perposition of propagating sine waves whose dispersion 
relation is close to the LDR (j3j) [8]. These observations 
rule out the possibility of strongly non linear structures 
being dominant in the motion of the plate as suggested 
by Boudaoud et al. [7 . 

The reason for this disagreement between theory and 
experiment remains an open question. Among the possi- 
bilities, one can list: 

1. finite size effects: these induce a quantization of the 
frequencies of the normal modes of the plate. If the 
spectral (frequency) width of the modes is not wide 
enough this could lead to discrete wave turbulence or 
laminated wave turbulence. The discretization of the 
frequencies may hinder the occurrence of 4-wave reso- 
nances and alter the energy flux [12 | [T3 l [T4]. 

2. deformation of the plate at rest: if the plate is not in- 
finitely flat at rest (which is always the case to some 
extend for real plates) then the non linear terms of the 
evolution equations may change from cubic to quadratic 
and this would also alter the energy transfer through 
3- wave resonances in place of 4-wave [7 . 

3. dissipation: if the dissipation is not localized at high 
frequency as assumed in the WTT, then it would al- 
ter the energy transfer. Energy would only be partially 
transferred among scales. Note that the source of dis- 
sipation here is not obvious and can involve the bulk 
material dissipation, radiation of acoustic waves, trans- 
fer of energy to the ambient air or to the support... It 
has been discussed to some extent in [7 . 

More recently, a Fourier Transform Profilometry (FTP) 
technique has been applied to this problem by Cobelli et 
al. [9 . This technique, coupled with a high speed camera, 
enables the measurement of the deformation of a signifi- 
cant portion of the surface of the plate, resolved in time 
and space. It was then possible to compute the space-time 
(k, uo) spectrum of the deformation. Here we report more 
extensive measurement and analysis using the FTP tech- 
nique with an upgraded camera. 



2 Description of the experiment and the data 
processing 

2.1 Experimental setup 

The experimental setup is similar to that described in [8j 
[9]. The plate is a stainless steel plate of size 2x1 m^ and 
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Fig. 1. Schematics of the experiment: the plate is 2 x 1 m^ 
and hangs under its own weight under a beam. The motion is 
forced by a electromagnetic shaker excited by a sine current at 
30 Hz which amplitude can be tuned. A high definition video 
projector projects a sine intensity pattern on the plate. The 
deformation of the plates induces a deformation of the pattern 
seen by a high speed camera. The field of view of the camera 
is about 0.8 x 1.1 m^. 



thickness 0.4 mm. It is bolted on a beam at the top along 
one short edge. It hangs freely under its own weight. At 
40 cm above the bottom, a electromagnetic shaker is an- 
chored on the plate and forces its motion (see schematics 
in fig. [1]). The excitation current is a pure sine at 30 Hz. 
Its amplitude can be tuned to explore various forcing in- 
tensities. The average power transmitted to the plate by 
the vibrator has been observed to be quadratic in the am- 
plitude of the excitation current [7 . Here we measure the 
average input power through arbitrary units of the forcing 
current in the interval [0.5,6] corresponding to a average 
power spanning an interval proportional to [0.25,36], i.e. 
two orders of magnitude. 

The FTP technique requires the projection of a sine 
intensity pattern /(x, y) cx sin(27rx/p) on the plate jl5[|16j. 
The plate is painted matte white to behave as a diffusive 
screen. Here x is the vertical coordinate. The wavelength 
of the pattern on the plate is p = 6.6 mm. The pattern 
is projected by a High Definition Epson video projector. 
The pattern is recorded by a Photron SAl high speed 
camera at a frame rate reaching 8000 frames per second 
(fps) at the highest forcing. The field of view of the camera 
is about 0.8 x 1.1 m^. This region is located about 10 cm 
away from the edges and in the central region between 
the top of the plate and the anchor point of the shaker 



in order to avoid possible boundary inhomogeneities. The 
field of view is significantly larger than reported in [9^ and 
the camera is both faster and more sensitive as that used 
previously. The FTP scheme is identical to that of [9l 
\Tb\. The projector and camera lenses have parallel optical 
axes. They are L = 2.17 m away from the plate at rest 
and the separation of the optical axes is D = 30 cm. 



2.2 Data processing 

The image processing is similar to that described in [15]. 
The deformation of the plate appears as a phase modula- 
tion of the pattern recorded by the camera. The phase of 
the deformed pattern is extracted and the deformation is 
computed as in [15^. In this way, one obtains a movie of 
the deformation of the plate over about one square me- 
ter. At each forcing amplitude, the recording rate of the 
camera is tuned to ensure that the Shannon criterium is 
satisfied up to frequencies at which the signal falls below 
the noise level. The velocity is computed by difference be- 
tween successive deformation frames. A snapshot of defor- 
mation and velocity is shown in fig.[2j The spectrum of the 
deformation being quite steep (about the deformed 

surface is quite smooth. On the contrary the spectrum of 
the velocity is quite flat (see below) resulting in a velocity 
that varies very quickly in space. 

From the deformation or velocity movies many analy- 
ses can be performed. As shown below, a Fourier transform 
in space can be applied to get the spectral wave ampli- 
tude v(k,t). This component is the sum of 'U+(k, and 
'U~(k, t) respectively the wave components with positive 
frequency cu > and the wave component with negative 
frequency (these two waves are traveling in opposite direc- 
tions). These two components can be separated numeri- 
cally. By symmetry one gets 'U+(k, t) = 'U~(— k, t)* (where 
stands for the complex conjugate) so that one only has 
to study (k, t) over the full complex plane to get all the 
information on waves traveling in all directions. Statistics 
of these complex amplitudes (distributions, correlations...) 
can be computed. For instance, the 2D space spectrum 
Ey{]i.) is simply the average (|'U+(k, t)p) where the aver- 
age is taken over time. By integrating the 2D spectrum 
over angles of k one gets the ID space spectrum Ey{k). 
A Fourier transform in time can be applied also, leading 
to the spectral wave amplitude v+(k, cj). The space-time 
Fourier spectrum £^^(k, cj) is obtained from the average 
(|'L'+(k, cj)!^) and it can also be integrated over the direc- 
tions of k to get Ey{k^uj). 



3 Spectral wave amplitudes 

3.1 Distribution of the wave amplitudes 

Figure [3] displays the measured distribution of the squared 
magnitude |'U+(k, t)p for values of k spanning the inter- 
val between the forcing wavenumber up to wavenumbers 
corresponding to a signal over noise ratio close to 1. The 
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(a) 




(b) 







X(m) 



Fig. 2. Example of the plate deformation and velocity 



distribution is exponential (Rayleigh distribution) as ex- 
pected for a gaussian complex amplitude v(k, t). No depar- 
ture is visible within the range of probabilities accessible 
with the current dataset, neither at small k nor at large 
values. Figure |4] shows the distribution of the phases of 
'L'(k, t) for the same set of values of k. The phases are cor- 
rected from the linear frequency, i.e. a linear drift a;L(k)t 
is subtracted. The distributions are flat at all values of k. 

The phase and the magnitude of v+(k, t) are also ex- 
pected to be independent variables. Figure [5] analyses the 
joint distribution of the magnitude and the phase for a 
given value of k. The left part of the figure shows the joint 
distribution P((/)(k, t), |'U+(k, (where (/)(k, t) stands for 
the phase of i;+(k, t)) and the right part shows the ratio 

p{<p{k,t),\v+{k,m 

P{<P{k,t))P{\v+{k,t)n- 

This ratio should be unity for independent phase and mag- 
nitudes. This is indeed observed in our experiment. Note 
that to test the hypothesis (v) on random phases of inter- 




v^(k)r/<|v^(k)| > 



Fig. 3. PDF of the squared magnitude of the spectral defor- 
mation components v~^(k,t). Here, k = (^,0) with values of 
kx equally spaced between 2.2 and 25. The various curves have 
been shifted vertically for clarity. Dashed lines correspond to 
exponential decay. 
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Fig. 4. PDF of the phase of the spectral deformation com- 
ponents v~^(k,t). Here, k = {kx,0) with values of kx equally 
spaced between 2.2 and 25. All curves are superimposed. 
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Fig. 5. example of the joint PDF of the magnitude and the 



phase of v^{k,t) for k/27r = (15,0). Left: P( 
ratio P(0, \v^\^)/P{(l))P{\v^\^) 



right: 
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acting modes, one would have to use higher order statis- 
tical tools such as bi- and tricoherence. 



3.2 Correlation of the wave amplitudes 



1/^1^=10.3 m"'' 



10 20 
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{v+(]^i,t)v+{]^2,tr) 



Fie. 6. Magnitude of the correlations . 

V(^+(kl,t)2)<^+(k2,t)2) 

Here, ki is chosen on the diagonal with its norm given in the 
title. 



Another test of the validity of the bases of the WTT 
is the fact that the spectral component should be uncor- 
related i.e. 



(^+(ki,t)^+(k2,tr) (xE,(ki)5(ki-k2) 



(6) 



This relation is tested in figure [6l The correlation are in- 
deed zero except on ki . The width of the correlation peak 
is due to the finite resolution of the space Fourier trans- 
form due to the size of the recorded picture. 



4 The space-time Fourier spectrum 
4.1 The Fourier spectra 

As the full space and time resolved deformation can be 
measured using the FTP technique, we gain access to var- 
ious Fourier spectra. Figure [7] shows the full space and 
time spectrum £^^(k, a;). It is very clear in this figure that 
energy is localized on a surface in the (k, cj) space. This 
is what is expected for waves with a dispersion relation. 
The cuts of EyCk^uj) at constant uj (fig.[8|) shows that en- 
ergy is localized on circles. The spectrum is shows some 
anisotropy at small uj due to the forcing that propagates 
waves preferentially towards the top of the plate. At higher 




ky/2pi[m ■] 



y2pi [m-'] 



Fig. 7. The space-time spectrum E(\i,uj) for P = 16. 3D view 
of the spectra with cuts at = 0, A;^ = and / = 500, 1000, 
1500, 2000 and 2500 Hz. Energy levels are color coded in log 
scale. 



frequencies, the energy is distributed continuously on the 
circles (except at the highest frequencies, see part [5] be- 
low) which shows that wave turbulence transfers energy 
in a continuous and almost isotropic manner in k space. 

Figure [9] displays the angle integrated spectrum E{k^ uj). 
The linear dispersion relation is shown as the black line. 
The line of maximum energy is localized very close to the 
linear dispersion relation (this point is analyzed further in 
the next section). 

The isotropy of the space spectrum is clearly seen in 
figure [TOTa). At low k the anisotropy of the forcing is obvi- 
ous but at higher k the lines of constant energy are circles. 
Angle integrated spectra E{k) are shown in fig. [TOTb) for 
P ranging from 0.25 to 36. The development of a short 
plateau (almost horizontal) can be seen that corresponds 
to an inert ial regime. This inert ial regime is not so wide 
so that hypothesis (iii) on scale separation of forcing and 
sinks can be questioned here. 

This inertial regime is much wider in time spectra 
(fig. [TT]) . due to the curvature of the dispersion relation 
UJ (X k'^. A f-^-^^ region is clearly visible as was already 
reported in previous work [71IH] • Using the linear dispersion 
relation this translates into k~^-^ as seen in fig. [ToT b). 



4.2 The non-linear dispersion relation 

A non linear dispersion relation (NLDR) can be obtained 
by extracting the position of the crest of the energy line in 
fig- El An example is shown in fig. [12] at the highest forcing 
P = 36. At high frequency it recovers the k^ scaling of the 
linear dispersion relation in the range of scales in which 
the spectra are decaying at the fastest rate. This allows 
actually to compute very precisely the constant C in the 
LDR ujl = Ck'^. Some uncertainty exists in the physical 
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Fig. 8. The space-time spectrum ^(k,C(;) for P — 16: cuts of 
the spectrum at given values of / (specified in the title of each 
subfigure) . 



parameters a and so on which provides estimations of 
C only within a few percent at best. Here we estimate 
C = 0.641 m^s~^ which gives the best agreement with 
the high k evolution of the NLDR. The accuracy on the 
estimation of C is about 0.2%. 

The WTT predicts a small correction to the LDR un- 
der the form ujNL{k) = ojL{k) + Auj{k) [2 . The shift 
Auj{k) is shown in fig. [TST a). It rises from zero and peaks 
at relatively low k (but still larger than the forcing at 
k = 2.S m~^. It decays to zero at the highest k avail- 
able (values of k at which the signal over noise ratio is 
larger than 1). The current setup with a fastest and more 
sensitive camera as in [9] enables to have a higher preci- 
sion of the measurement of Auj and C which must be the 
reason for the difference with the data shown in [9 . The 
curves at various P are self similar and collapse to a master 
curve when plotted Acu/P^^'^ versus k/P^^^ (fig- [Hb)). 
The scaling of k/P^^^ corresponds through the LDR to 




00 -1500 1500 3000 
frequency (Hz) 

Fig. 9. The angle integrated space-time spectrum E{k,(jj) for 
P = 16. 



the scaling f /P^^^ that is also shown to collapse the time 
spectra [3[8]. It is due to the P^/^ scaling of the cutoff 
frequency observed in the time spectra (fig. (TTJ [8 ). The 
physical origin of the cutoff frequency is still unclear but 
it could be independent of dissipation as such a cutoff fre- 
quency can be exhibited by dimensional analysis of the 
FAUppl-Von Karman equations ([I])&(|2]) [8 . The analyt- 
ical expression of Auj has not been derived yet from the 
WTT in the case of the elastic plate. Even if the WTT 
does not agree with the observed spectra, the scaling ob- 
served for Auj may still be related to theoretical predic- 
tions. If one considers Auj as a measure of the degree of 
non linearity, it is the strongest slightly above the forcing 
and decays to zero at high frequencies. Thus we do not 
expect the hypotheses of weak non linearity to breakdown 
at high frequencies as is the case in some systems |2] . The 
shift to the LDR is about 50% of the LDR at the forcing 
scales and much smaller in the inert ial range. Together 
with the absence of harmonics, this observation can be 
considered as a indirect validation of hypothesis (ii) of 
weak non linearity. 



4.3 The width of the dispersion relation 

The theoretical correction Auj can be complex a priori. 
Thus one expects that the width of the energy line around 
the NLDR in figure [9] should be increased by non linear 
effects due to the imaginary part of Auj . This width has 
been computed in two ways in fig. [Ml either a /c- width at 
constant / or a /-width at constant k. The /c- width has 
a lower bound due to the resolution of the space Fourier 
transform on the size of the picture. It is shown as a dashed 
line in fig. [TU The measured /c- width is indeed larger than 
this lower bound and is seen to increase monotonically 
with P. The lower bound of the /c-width can be translated 
into a lower bound for the /-width through the group 
velocity of the LDR. It gives the dashed line in the insert 
of fig. [m which is indeed a lower bound for the measured 
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(a) 



k I2n (m"'') 




(b) 
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Fig. 10. The space spectra ^(k) (a) for P ■ 
integrated space spectrum E{k) (b) for P ■ 



16 and the angle 
0.25, 1, 4, 9, 16, 



25 and 36 (from bottom to top). The dashed hne is /c " . The 
forcing is at /c = 2.8 . 



/-width. The width behaves similarly to the shift to the 
LDR, with a peak slightly above the forcing wavenumber 
and a decay at the largest frequencies. 



5 Finite size effects 

The hypothesis (i) of WTT assumes an infinite system 
to avoid discrete values of the frequencies of the available 
vibration modes. Figure [15] shows a cut of the spectrum 
^(k, cj) in the plane kx = 0. At low frequencies, the en- 
ergy line is continuous but at the highest frequencies the 
energy is localized in peaks. These peaks appear in the re- 
gion of the frequency spectrum in which the spectrum is 
decaying exponentially and is the weakest. In this region 
the discreteness of the normal modes of the plate seem to 
appear. They are the most visible in the vicinity of the 
plane /c^^ = in the (k, uo) space. The y direction is along 
the short side of the plate. The separation of the modes is 



LU 




Fig. 11. The time spectrum E{uj) for P = 0.25, 1, 4, 9, 16, 
25 and 36 (from bottom to top). The spurious peaks in the 
spectra are either harmonics of the forcing or due to flickering 
of the videoprojector acting at small k as also seen in fig. [9] 
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Fig. 12. The non linear dispersion relation at P — 36. The 
insert is in loglog scale. 



the widest along this direction. One can extract the am- 
plitude of the energy line along the NLDR of fig. [151 The 
result is shown in fig.[T6]for various forcing amplitudes and 
in the frequency window / G [500,900] Hz. The spectral 
width of the peak is increasing with the forcing until the 
peaks merge at the highest forcing. In the same frequency 
window, similar peaks can be seen also in other directions 
of k {ky = for instance) but the interval between peaks 
is much lower so that the peaks merge for much smaller 
forcing amplitudes. This can be observed in the cuts of 
£^(k, cj) for given values of the frequency shown in fig. [51 
The two lower cuts show clear localized peaks of energy 
near kx = ^ (that are not due to insufficient convergence 
of the average). When kx increases the energy is continu- 
ously distributed along the circle that corresponds to the 
dispersion relation. 
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Fig. 13. (a) The shift of non hnear dispersion relation in re- 
spect to the hnear dispersion relation Auj{k) — uJNL{k)—ujL{k) 
for P = 0.25, 1, 4, 9, 16, 25 and 36 (from bottom to top). The 
forcing is at /c = 2.8 m"^ (b) Au{k)/P^^'^ versus k/P^^^. 



Finite size effects are visible only in the high frequency, 
fast decaying region of the spectra. In this region, the mea- 
surement of the shift to the linear dispersion relation show 
that the non linearity is the weakest. The widening of the 
discrete frequencies is thus not strong enough to observe 
a continuous spectrum. In the power law region of the 
spectra, the energy is distributed continuously as the non 
linearity is stronger. 

In [17 , it was suggested that finite size effects could 
cause a "sandpile behaviour" of the energy flux that even- 
tually leads to a different scaling of the spectrum. The 
argument is the following: energy will be able to flow 
across scales only when the amplitude of the waves in- 
duces a widening of the frequencies of adjacent (in k- 
space) modes equal to the separation of their frequencies. 
The frequency widening of the mode is of order and 



it should be equal to the mode separation which is of or- 
der ~ k/L (with L the size of the finite container 
which is of order 1 m here). From [6 the collision integral 
scales as uj{kY~^^ for a spectrum n{k) ex uo{k)~^ . The 
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Fig. 14. The A;-width of the energy line of fig. [9] (see text) 
for P = 1, 4, 9, 16, 25 and 36 (from bottom to top). The 
dashed line is the resolution of the space Fourier transform. The 
insert shows the /-width. The dashed line is the resolution of 
the space Fourier transform translated in frequency space (see 
text). 
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Fig. 15. Cut of the space time spectrum £^(k,cj) at kx — 0. 
Discrete frequencies are clearly visible. 



so-called "critical spectrum" is then obtained by solving 
uj(kY~'^^ (X cj(/c)^/^ which provides x = 1/4. In [6 , the 
collision integral is said to converge for 0.5 < x < 1.2. Here 
we get X = 1/4: which does not ensure the locality of the 
interactions in Fourier space (hyp. (iv)). The Kolmogorov- 
Zakharov solution corresponds to x = 1 i.e. steeper than 
the possible critical spectrum. The observed spectrum is 
even steeper than the Kolmogorov-Zakharov prediction 
and thus not in agreement with the critical spectrum pre- 
diction. As the latter is affected by the non locality is- 
sue, this deserves more theoretical investigations. What is 
observed in fig. [15] is that the non linear widening of the 



Nicolas Mordant: Fourier analysis of wave turbulence in a thin elastic plate 



9 



Q 
(/) 
Q_ 



10' 



10' 



10' 




500 



600 



700 
f[Hz] 



800 900 



Fig. 16. Evolution of the spectral energy on the dispersion 
relation in the plane kx = for P = 1, 4, 9, 16 and 36 (from 
bottom to top). The curves have been shifted vertically by a 
factor 1.5 for clarity. 



modes is greater than the mode separation at low frequen- 
cies so that a Kolmogorov-Zakharov cascade may proceed. 
When the widening is of the same order as the separation, 
avalanches may occur. In the highest part of the spectrum, 
the modes are apparent and it must correspond to another 
phenomenology. 



6 Conclusions 

All measurements of wave turbulence in an elastic plate 
reported in this article are qualitatively consistent with 
the kinetic Weak Turbulence Theory. A transition from 
discrete wave turbulence to a kinetic wave turbulence is 
observed when the non linearity is sufficiently strong, as 
described by Kartashova [12,13,14 , with a intermediate 
mesoscopic or laminated turbulence. The existence of a 
frequency cutoff leading to an exponential decay of the 
spectrum at large frequency may be due to this transition 
to discrete turbulence. It has been observed that the spec- 
tra are steeper in mesoscopic or discrete turbulence due 
to a limited set of available resonant frequencies. 

The non linearity is weak in this system as quantified 
by the small shift of the observed non linear dispersion 
relation in respect to the LDR. As expected from the phe- 
nomenology, isotropy of the energy spectra is restored at 
large frequencies despite the anisotropic forcing. Despite 
this agreement with the hypotheses of the WTT, the en- 
ergy spectrum is not in quantitative agreement with the 
theory. Recent work by Connaughton [5] have shown that 
in numerical simulation, the truncation scheme can induce 
a positive or a negative bottleneck effect that can be ex- 
tended in a quite large interval of frequencies. The dissipa- 
tive mechanism are complex and unknown to some extend 
in the case of the elastic plates. Nevertheless a negative 
bottleneck effect could be responsible for the quantitative 



disagreement observed here especially since the interval 
of accessible frequencies (let alone wavenumbers) is not 
extremely wide. 

The author thanks the french "Agence Nationale de 
la Recherche" for its funding under grant TURBONDE 
BLAN07-3-197846. 
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